Ftactional-flux vortices and spin superfluidity in triplet superconductors 
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We discuss a novel type of fractional flux vortices along with integer flux vortices in Kosterlitz- 
Thouless transitions in a triplet superconductor. We show that under certain conditions a spin- 
triplet superconductor should exhibit a novel state of spin superfluidity without superconductivity. 



Superconductors with triplet pairing allow for a rich 
variety of topological defects and phase transitions P, 0, 
S 0) S 's'l . For example in p-wave superconductors there 
exist different realizations of fractional vortices such as: 
the half-quantum vortex (the Alice string) , which is a vor- 
tex where phase changes by tt and spin is reversed when 
we go around the vortex core, fractional vortices trapped 
on a grain boundary, fractional flux trapped by twisted 
wire and other possibilities (for an excellent review see 
Another subject which was intensively studied is 
the split phase transitions in triplet systems (e^. in the 
presence of disorder or magnetic field) '5, '3, There 
have also been particularly interesting studies of split 
transitions in a neutral p-wave superfluid in connection 
with thin films of liquid ^He 0i @|- Also the related 
questions of various partial symmetry breakdowns are 
relevant and were studied in spin-1 Bose condensates in 
optical traps 0| and bilayer Quantum Hall systems '^l- 
In this Letter we discuss the effect of flux fractionaliza- 
tion due to spin-orbit coupling, as well as the appearance 
of neutral vortices in triplet superconductors and discuss 
its influence on the Kosterlitz-Thouless (KT) transitions. 
Based on topological arguments we predict the existence 
of a novel type of ordering in spin-triplet system - the 
spin-superfluid nonsuperconducting state. 

We consider a model spin triplet superconductor simi- 
lar to one considered in with order parameter ^'a(r) — 
■\/n(x)Ca(x) where (a = 1,0,-1) and C is a normalized 
spinor C,^ ■ Q = I. If one neglects mixed gradient terms, 
the free energy density can be written as 0, 0, 0| : 
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where < F QFab(b is the spin. The pres- 

ence of a crystal lattice makes the mass M a sym- 
metric tensor Mij. Spinors are related to one an- 
other by gauge transformation e*^ and spin rotations 
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Euler angles. The matrices {F^, Fy, F^} are 
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The ground state structure of 5'a(r) can be foun d by 
minimizing the energy with fixed particle number |5|, [llj . 
When C2 < the energy is minimized by < F >^= 1 and 
the ground state spinor and density are [ll| 
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In this case there exists an equivalence between the gauge 
transformation associated with 6 and spin rotations asso- 
ciated with r and the vacuum manifold is S'0(3) [l|l[lll. 
The equation for the supercurrent in this case depends on 
phase gradients and also on spin texture (the Mermin-Ho 
relation) : 
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In it was discussed that can be rewritten 

in the form of an extended Faddeev model re- 
lated to that derived for two-gap superconductors 0, 
[l6i |. That is, introducing notations — -£r, s = 

(sin /? cos a, sin /3 sin a, cos /3), C = and separating 

variables, the free energy density can be rewritten 
in terms of only gauge-invariant variables as 
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Let us now allow for a spin-orbit coupling breaking the 
0(3) symmetry. In the simplest case one should add to 
© a term 



(6) 



The case when there is only one nonzero coefficient 
7zz < corresponds to the "easy-axis" situation, then 
0(3) spin symmetry is broken down to a point. When 
> we have the "easy plane" case and the symmetry 
is broken down to a circle. Adding the terms "fxx = Jyy 
still leads to an easy-plane or an easy-axis case. However, 
in the general case there should be higher order correc- 
tions e.g. fourth order in Si. Such corrections shift the 
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ground state from easy-axis or easy plane and the ener- 
getically preferred state of spin direction can be a circle 
on the unit sphere as shown on Fig. 1. Such a situ- 
ation is known in many magnetic systems |l7l |. Let us 
now examine the consequences of the spin-orbit coupling 
for the vortices of map. Let us consider the 

case when spin-orbit coupling breaks spin rotation sym- 
metry down to an arbitrary circle on 5^ characterized by 
f} = (3o (see Fig. 1) (we consider the range < /3o < tt). 
One can observe that when (3 — (3q — const we have 



FIG. 1: Ground state of s for Pq / 7r/2 

s • ViS X VjS cx sin/3[Vi/3VjQ; — V jfiVia] — 0. In this 
case, in the London limit, the eq. I© becomes: 
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From lO it follows that the system allows the following 
vortices of map: (i) A vortex where 9 (or r) 

changes by 27r around the core. This is the analogue 
of Abrikosov vortex in ordinary superconductors and it 
carries one flux quantum, (ii) A vortex where a changes 
by 27r around the core. This vortex has a neutral vorticity 
with the stiffness |^ sin^ (3o described by the first term 
in |7J) and also there is a supercurrent around this vortex 
due to Mermin-Ho relation, which is given by 
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From here it follows that the magnetic flux of this vor- 
tex is determined solely by spin-orbit coupling and is an 
arbitrary fraction of magnetic flux quantum: 
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where cr is a path around the core at a distance much 
larger than A and $o is the magnetic flux quantum. A 
special situation appears in the "easy-plane" case, when 
/3o = 7r/2, and a vortex in spin degree of freedom does 
not carry magnetic flux. 

Let us now discuss the consequences of the above 
topological excitations to phase diagram of a quasi-two- 
dimensional model spin-triplet superconductor where the 
0(3) symmetry is broken to 0(2). We consider the case 



when this symmetry is broken strongly and restrict dis- 
cussion to the range < /3 < 7r/2, with the main focus 
on the case /3 — tt/2. Thus for simplicity we do not con- 
sider in this Letter the case of two degenerate minima 
with f3 — Po and /3 = tt — /3o). From now on we discuss 
a quasi-two dimensional system of thickness d where the 
effective magnetic field penetration length A = Ap is re- 
lated to Ginzbur g-La ndau penetration length Xql via the 
Pearl expression [ijj Ap = XQi^/d. 

We begin with the simplest limit when the magnetic 
field penetration length A is larger than the sample size 
and Po — 7r/2 (the easy plane case). In such a sit- 
uation we can neglect coupling to the vector poten- 
tial in iQ. Then we can observe that vortices where 
{A9 = <f dlVO = 27r,Aa = § dlVa = 0) are mapped onto 
vortices in the ordinary one component U{\) model with 
stiffness Jg = while the vortices (Aa = <fdlVa = 
2TT,Ad — ^ dlWd = 0) are mapped onto vortices in a 
C/(l)-model with stiffness Jq — Thus in the limit 

A — !■ CX) there are no corrections to vortex interaction due 
to the Meissner effect, and the system undergoes two KT 
phase transitions at the temperatures: 
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(equations for Tkt should be solved self-consistenly with 
the equations for gap modulus, see e.g. the review 0|). 
This splitting of the phase transitions corresponds to sep- 
arate onset of quasi-long-range order in phases 6 and a. 

The regime of short penetration length is however prin- 
cipally different and more interesting. In this limit the in- 
teraction which is mediated by charged current between 
vortices and antivortices {A6 — ±27r) is exponentially 
screened at the length scale A and such vortices do not 
undergo a KT transition. Thus the variable 9 is disor- 
dered and the system is not superconducting. However 
in this limit the vortices with (Aa = ±27r) still have a 
long range interaction mediated by the neutral mode and 
characterized by stiffness: 
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(as follows from the first term in {T))). Thus the "spin 
vortices" in the regime 7r/3 < Pq < 7r/2 can undergo a 
true KT transition at the temperature 
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this phase transition is associated with establishment of 
quasi-long range order in the phase variable a. Note that 
in this state the superconducting phase 9 is still disor- 
dered so is not a superconducting phase transition 
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but is a transition to a "spin superfluid state^\ Let us 
emphasise that the vortices with (Aa = ±27r) have non- 
trivial contribution in the second term in 10 and also 
carry a fraction of flux quantum @. 

We should observe that the situation when < /3o < 
7r/3 is very special. That is, in the short- A limit the long 
range interaction of vortices is mediated by the neutral 
current (coming from the neutral mode described by the 
first term in Q ) and besides that there is a finite energy 
contribution from the charged current (described by the 
the second term in {T))). Although the later contribu- 
tion is irrelevant for vortex interaction (in Coulomb gas 
mapping it plays the same role as the core energy), how- 
ever the charged mode is important in the following re- 
spect: In the regime < Pq < n/3 the composite vortices 
(Aa — ±27r, A^ = ±27r) have smaller energy than ele- 
mentary vortices (Aa = ±27r, A9 = 0) since the former 
vortices carry a smaller fraction of magnetic flux quan- 
tum and thus have smaller kinetic energy of charged su- 
percurrent. Thus when < /3o < tt/S there will be a KT 
transition of the composite vortices (Aa = 2tt, A9 = 2tt) 
and antivortices (Aa = —2%, A6 — —2tt) at the tem- 
perature H12() . This is a rather unique feature, rooted 
in the nontrivial coupling of phases by vector poten- 
tial in the triplet superconductor, when a KT transi- 
tion is governed by topological defects with high topo- 
logical charge. Note that in this regime the system al- 
lows also "purely charged" integer flux vortices (that is, 
without neutral superfiow) (Aa — 0, AO = ±27r). These 
vortices have only screened interaction. In this regime 
they are always liberated and disorder superconducting 
phase. Thus, also in the regime < /?o < tt/S the phase 
transition H12|l is a transition to a nonsuperconducting 
spin-superfluid state in spite it is mediated by composite 
vortices which have topological windings in both phases 
(Aa = ±27r, A6I = ±2tt). 

Let us now discuss qualitatively regimes of intermedi- 
ate penetration length. In a general case of large but fi- 
nite A, the interaction of vortices and antivortices, which 
is mediated by charged current, receives correction due to 
Meissner screening (for a detailed review and citations see 
^3). Then, if A is sufficiently large, a gauged purely U{1) 
system, althought, strictly speaking, does not display a 
KT transition into a superconducting state, is however 
believed to undergo a "would be" KT crossover 0. So, 
Tkt at finite A is transformed into a characteristic tem- 
perature of a broader crossover and also is shifted down 
comparing to the situation when there is no Meissner 
screening 19]. In the system {TJ, at finite- A, the vortices 
{AO — ±.2tt) interact via charged mode with the effective 
stiffness: 
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and with charged mode, which effective stiffness in the 
case 7r/3 < /3o < 7r/2 is 

^A.V3<;3o<V2(r) = eos^ /3o^^ (15) 

When < /9o < 7r/3 the relevant topological excitation 
is the composite vortex (Aa — 2n, AO = 2tt) for which 

^A,0</5e<V3(r) ^ [1 _ COS/3o]^^^ (16) 

Let us now consider the easy-plane case Pq = tt/2. 
In type-I limit, as discussed above, there is only spin- 
superfiuid KT phase transition while superconductivity 
sets in only at T ^ 0. When we increase penetration 
length the system undergoes a washed out superconduct- 
ing KT crossover [lij at some nonzero temperature. In 
the extreme limit A — > oo we arrive at the above limit of 
two true KT transitions and an intermediate state where 
there is no quasi-long range order in spin degrees of free- 
dom while the system is superconducting. The case of ar- 
bitrary f3a and A appears to be a very interesting subject 
for numerical simulations. A schematic phase diagram in 
the simplest case Pq = tt/2 is given on Fig. El 
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while the vortices with (Aa = ±27r) interact by means 
of both the neutral spin mode with the effective stiffness 



FIG. 2: A schematic phase diagram of the system in the 
simplest regime /3o = it/2. The temperature of onset of su- 
perconductivity, at finite penetration length is a characteristic 
temperature of a washed out crossover. 

Let us discuss the physical interpretation of the "spin- 
superfiuid" state. In this state a quasi-long-range or- 
der is retained in the spin angle a. We observe that 
all three complex components in Eq. Q have a pref- 
actor e**. Thus in the spin-superfiuid state, all three 
phases of the spin components are disordered. Con- 
sequently, there is no dissipationless supercurrent 10} . 
However the system retains hidden order, that is, taking 
the difference of the total phases of the first and third 
components, the contributions from the phase is can- 
celled, while the contribution of the ordered phase a is 
not. That means that there is a quasi-long-range or- 
der in a composite order parameter. Such a composite 
order parameter can be identified though the principle 
that a squared modulus of its gradient should give the 
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kinetic term of the composite neutral mode in the GL 
functional after separation of variables (UJ. This term 
is the following combination of the fields (a which is de- 
coupled from the vector potential but is gauge invariant: 

(?i2n/2M)[VCaVC + (l/4)(C:VCa - CaVC)2]|/3=/5o This 
corresponds to the following order parameter for the com- 
posite neutral mode: 



1 



sin /3oe*' 



(17) 



In the hydrodynamic limit the spin superfluid velocity 
is given by V5 = sin^ /3oVa. This expression has 
a physical interpretation: in the easy-axis limits (when 
/3o = or /3o = tt) there is no neutral U{1) symmetry 
and we find that also J5 — 0. On the other hand in the 
easy-plane case when Pq — ^ we have l^"!] = l^*-!! = 
|\l/o|/-\/2. We note that is neutral with respect to a 
and thus does not participate in spin superfluidity. This 
gives a natural explanation for the extra factor 1/2 in 
front of the phase gradient in the expression for Vg in the 
easy-plane situation, compared to the similar coeffcient 
in front of W0 in the expression for (l/2e)J given by eq. 

O- 

The configuration of the phases of the three complex 
components of the order parameter in the spin-superfluid 
state is illustrated in Fig. 13 

SPIN-SUPERFLUID STATE 




FIG. 3: A schematic illustration of the configuration of 
phases of individual spin components in different points in 
space in the spin-superfluid state. 



transition by a FN measurement, in subsequent measure- 
ments of the conductivity or application of an external 
field, the system will not behave as a superconductor {in) 
If the corrections shifting the spin from the equator on S*^ 
are absent - there nonetheless is a possibility to observe 
the transition to the spin-superfluid state if an external 
field is applied. Then the Zeeman effect would result in 
a term linear in Sz in @ thus shifting spins from equator 
on , consequently the spin vortices will acquire a frac- 
tion of flux quantum and the transition will be detectable 
via FN measurement. 

In conclusion, a triplet superconductor is a multicom- 
ponent charged system, where the vector potential is cou- 
pled to phase variables in a nontrivial way which results 
in the existence of neutral modes. The nontrivial cou- 
pling to vector potential of the components of the or- 
der parameter becomes particularly important in (quasi)- 
two dimensions where the system undergoes topological 
phase transitions. We considered a simple model of a 
spin-1 superconductor and have shown that topological 
considerations lead to a novel physical state: the "spin- 
superfluidity" when only an opposite superflow of two of 
three spin components is allowed. The results hold true 
in the simplest case of the easy-plan situation, however 
we extended the discussion to a general case when the 
ground state is shifted from the equator on 5^ which, as 
we have shown, leads to flux fractionalization. We also 
pointed out that there is a range of parameters when 
the quasi-long range order sets in only in the supercon- 
ducting phase and pointed out an unusual KT transition 
mediated by defects of high topological charge. 
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Possible experimental probes of the spin-superfluid 
state are: (i) When there are corrections shifting spin 
from the easy plane, the vortices in the a-phase carry a 
fraction of flux quantum. The intrinsic flux noise of a dc 
SQUID is about 10^^ ^q/\/Hz which makes flux-noise 
(FN) measurements sensitive to spin-superfluid transi- 
tion even if a tiniest "shifting correction" is present. An 
observation of a KT transition in the spin-superfluid state 
in an FN experiment should allow to distinguish it from 
a superconducting transition: Because of the existence 
of the neutral mode, the transition should be sharp and 
a true jump in superfluid density should be recoverable 
from flux-noise data, in contrast to similar measurements 
of wider KT crossovers in ordinary superconductors, (ii) 
Upon the transition to spin superfluid state a "true" 
superconductivity is absent. So after detection of the 
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